A deep submillimetre survey of the Galactic Centre 



D. Pierce-Price^, J. S. Richer^, J. S. Greaves^, W. S. Holland^, T. Jenness^, A. N. Lasenby^, 
G. J. White^'^ H. E. Matthews'''^, D. Ward-Thompson^, W. R. F. Dent^ R. Zylka^, P. Mezger^, 
T. Hasegawa^, T. Oka^^, A. Omont^^ and G. Gilmore^^ 

ABSTRACT 

We present first results from a submillimetre continuum survey of the Galactic Cen- 
tre 'Central Molecular Zone' (CMZ), made with SCUBA on the James Clerk Maxwell 
Telescope. SCUBA's scan-map mode has allowed us to make extremely wide-field maps 
of thermal dust emission with unprecedented speed and sensitivity. We also discuss 
some issues related to the elimination of artefacts in scan-map data. Our simultaneous 
850/450 //m maps have a total size of approximately 2.8 x 0.5° (400 x 75 pc) elongated 
along the galactic plane. It covers the Sgr A region — including Sgr A*, the circumnuclear 
disc, and the +20km/s and +50km/s clouds; the area around the Pistol; Sgr B2 — the 
brightest feature on the map; and at its Galactic Western and Eastern edges the Sgr C 
and Sgr D regions. There are many striking features such as filaments and shell-like 
structures, as well as point sources such as Sgr A* itself. The total mass in the Central 
Molecular Zone is greater than that revealed in previous optically-thin molecular line 
maps by a factor of ~ 3, and new details are revealed on scales down to 0.33 pc across 
this 400pc-wide region. 
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Subject headings: Galaxy: center — ISM: clouds — ISM: structure — dust, extinction 
— techniques: image processing — submillimeter 

1. Introduction 

The Galactic Centre provides our best opportunity to study the astrophysical processes in a 
galactic nucleus with high spatial resolution. These include star-formation and the molecular cloud 
structures in a region subjected to strong shear forces, magnetic fields, and gravitational potentials. 
Thermal dust continuum emission traces the temperature-weighted column density of material in 
a less biased way than molecular line maps, which are affected by excitation, optical depth, and 
abundance variations. Thus, we can use the continuum emission to map the mass distribution in 
molecular clouds with good accuracy and sensitivity. 

Wc used SCUBA (Holland et al. 1999) on the 15-m James Clerk Maxwell Telescope (JCMT), 
between April 1998 and April 2000, to map the Galactic Centre region in the submillimetre con- 
tinuum at 450 /im and 850 /xm. The new images cover a region 4 times larger and ~ 100 times 
deeper than earlier surveys (Lis & Carlstrom 1994). The map size is approximately 2.8 x 0.5°, 
or 400 X 75 pc at an assumed distance of 8.5 kpc. The map covers the entire 'Central Molecular 
Zone' (CMZ), which extends to radii of 200 pc and contains up to 10% of the Galaxy's molecular 
interstellar medium (ISM) (Morris & Serabyn 1996). 

2. Observations and data reduction 

SCUBA (Holland et al. 1999) is a submillimetre continuum bolometer camera with two arrays, 
one of 91 pixels optimised for 450 /nm, and one of 37 pixels optimised for 850 jim. In use, both arrays 
observe simultaneously. The resolution is 8" FWHM at 450 /im and 15" at 850 /im. 

SCUBA was used in scan-map mode, in which the array is scanned across the sky at a rate of 
24" per second whilst chopping the secondary mirror at 7.8 Hz. This chopping produces images in 
which the sky is convolved with the primary beam of the telescope and a dual-beam function, of 
positive and negative delta functions separated by the chop throw. 

The majority of the map was observed six separate times, with chop throws of 20, 30, and 
65" in both R.A. and Dec, although some fields at the Galactic Eastern and Western ends were 
observed with 30, 44, and 68" chops. Each dual-beam image contains no information at spatial 
frequencies equal to the inverse chop throw and its harmonics, and the constant power level (zero 
spatial frequency) is also lost. 

The data were reduced with SURF, the standard SCUBA reduction software (Jenness Sz Light- 
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foot 1999), using the "Emerson 2" deconvolution algorithm. This algorithm combines the multiple 
chop-throw images of the map fields. A weighted average of these observations is formed in Fourier 
space, and the reverse transform yields the final map (Jenness &: Lightfoot 1998). Although the 
combination of multiple chops improves the spatial frequency coverage, it is unavoidable that spa- 
tial scales more extended than a few times the maximum chop-throw are measured with reduced 
sensitivity, and the zero spatial frequency is still lost. Ultimately, the survey is sensitive to spa- 
tial scales in the approximate range 0.3 pc (the beam size at 450 /um) to ~ 10 pc (a few times the 
maximum chop-throw). 

Calibration was performed from maps of Uranus and Mars, and all data were corrected for 
atmospheric extinction using skydips. Zenith opacities were in the range 0.1-0.45 at 850 /Ltm and 
0.33-2.5 at 450 /im. Most of the data in the central 1.5° region was observed in extremely good 
weather, with r(850/im) < 0.13, r(450/im) < 0.5. 

To calibrate fluxes of extended structures, the integrated signal in a 1' diameter aperture 
was used to calculate a conversion factor in Jy arcsec"^ V~^. Calibration relative to the primary 
calibrators is estimated to be accurate to 5% at 850 /um and 20% at 450 /nm, based on the dispersion 
of measured flux conversion factors, although absolute cloud fluxes also depend on corrections for 
constant map offsets as described in section 2.1. 

The survey has resolutions of approximately 8" at 450 /um and 15" at 850 /xm, corresponding 
to distances of 0.33 pc and 0.62 pc at a distance of 8.5 kpc. The 1-cr sensitivities per beam are 
approximately 30 mJy and 300 mJy at 850 /xm and 450 fim respectively, or ~ 20 Mq per beam with 
the assumptions in section 3.5. 

2.1. Corrections for artefacts 

Making large-scale SCUBA scan-maps is difficult for several reasons. The total power level is 
not recorded, so there may be a constant offset on each field, and this offset may in general vary 
between fields in the mosaic. Within individual fields, one must account for baseline removal on the 
time-series signals from the bolometers. The dual-beam chop will suppress information at certain 
spatial frequencies, as described above, and as there is so much extended structure in the CMZ it 
is difficult to scan or chop onto regions of no emission. 

Some previous SCUBA scan-maps have sometimes exhibited reduction artefacts, such as neg- 
ative troughs around bright sources. We have found that negative troughs can be greatly reduced 
by careful baseline removal on a per-bolometer basis. A single 10' x 10' map is contructed from 
typically 10-15 scans of the array across the field. Calculating a baseline separately for each scan 
will produce large spurious offsets for those scans which pass over bright structure, so ideally it 
is best to determine a baseline using only those scans which contain little bright emission. This 
proved impossible for a region as rich in structure as the Galactic Centre, but we have found that 
using the median level over all the scans for each bolometer does produce good results without the 
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spurious offsets. 

Tfiis is tfie largest area ever scan-mapped at 450 and 850 /xm, witfi more than fifty overlapping 
10' X 10' fields, which must be mosaiced to form the final image. The Starlink CCDPACK package 
(Draper, Taylor, &: Allan 2000) was used to combine the dual-beam images and remove relative 
level offsets between fields of order 10 and 0.5 mJy arcsec"^ at 450 and 850 ;Um respectively. We 
then applied further corrections to ensure that the pixel mean calculated along the chop throw 
direction was zero, for all of the mosaiced dual-beam maps. This must be true, as long as the 
edges of the map are essentially emission-free, and is necessary because incorrect zero levels in 
the mosaiced dual-beam maps will translate to slopes in the final images. We finally corrected for 
the constant zero offset on the entire map by measuring the typical zero levels from map areas 
devoid of identifiable clouds, and assigning an uncertainty based on the level variations. The final 
uncertainty in baseline levels is 15 and 1 mJy arcsec"^ at 450 and 850 pLui respectively, the former 
being higher in part due to the greater calibration uncertainties. These uncertainties have been 
neglected, compared to the sensitivity noise levels of the maps. 

3. Analysis 

3.1. General features 

Our maps (figures 1 and 2) show good correspondence with previous surveys of the region, 
such as that by Lis & Carlstrom (1994) but with much higher sensitivity and resolution. At 850 //m 
we can clearly see the central source Sgr A* and its Circum-Nuclear Disc (CND) (/ = —0.05, 
b = —0.05)), although these features are less pronounced at 450 /xm. We also see the compact 
GMCs M-0.13-0.08 (the +20km/s cloud) and M-0.02-0.07 (the +50km/s cloud). Close to this is 
the high- velocity molecular cloud CO 0.02-0.02 (Oka et al. 1999). At higher galactic longitude we 
see the Pistol region {I = 0.15, b = -0.05), other GMCs such as MO. 25+0.01 (Lis et al. 1994), 
Sgr Bl = 0.5) and B2 (/ = 0.66), and the Sgr T) (I = 1.1) star-forming region. At more negative 
longitudes, the map extends to the Sgr C (Z = —0.5) star forming region. All of the field is rich 
in extended structure, dominated by filamentary clouds and cavities which may be associated with 
supernova remnants. Such structures have also been observed by Oka et al. (1998) in CO (J = 1-0) 
emission. 

3.2. Sagittarius A* 

At both wavelengths we have a clear detection of the Sgr A* point-source. By taking slices 
through the position of Sgr A* to determine and remove the local background, we estimate the 
flux from Sgr A* itself, with quoted 1-a uncertainties. The results are also subject to the overall 
calibration uncertainties described in section 2. At 850 /xm we measure a flux of 2.6 ± 0.3 Jy. This 
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is consistent with the value measured from CSO-JCMT interferometry (Serabyn et al. 1997) of 
3.2 lb 0.7 Jy. At 450 iim, the same procedure gives a flux for Sgr A* of 1.2 it 0.4 Jy. This is consistent 
with the upper limit of 1.5 Jy found by Dent et al. (1993). These fluxes are somewhat lower than 
those measured by Zylka et al. (1995), which were 3.5 ± 0.5 Jy at 800 //m (rather than 850 //m) 
and 3.0 it 1.0 Jy at 450 //m, suggesting that there may be some variability of the Sgr A* point 
source. For example, Tsuboi, Miyazaki, & Tsutsumi (1999) observed a flare at 3mm during which 
the flux doubled, and Serabyn et al. (1997) found a flux of 7 ± 2 Jy at 350 //m. Our fluxes indicate 
a spectral index for the emission from Sgr A* of — 1.5]'l]^ 4, consistent with a non-dust component 
in the emission. Flux values for the synchrotron component from the postulated black hole are 
discussed further by Aitken et al. (2000). 



3.3. Greybody fitting 

Although we can calculate spectral indices in the submillimetre regime from the SCUBA data, 
these data do not constrain the overall form of the dust emission, and do not allow us to estimate 
the temperature of the dust. To do this, we need data around the peak of the emission, which for 
dust at these temperatures is in the far infra-red. 

We have combined our submillimetre data with 100 jum data from the IRAS Galaxy At- 
las (Cao et al. 1997), and fitted a single temperature greybody spectrum of the form Ii, = 
(1 — cxp(— r))_B^(T), where the optical depth r is expressed in terms of an index (3 and a crit- 
ical frequency fcrit at which the emission becomes optically thick: r = (i^/t'crit)'^- Although the 
IRAS Galaxy Atlas also contains 60 /im data, we believe that emission at this wavelength traces a 
higher temperature component. By fitting to the 850 /xm, 450 //m, and 100 //m points we calculate 
properties of the cold dust only. 

We have performed a least-squares greybody fit to the three datasets over the survey area, 
fitting for temperature T and critical frequency fcrit- We fixed /3 = 2, the theoretical maximum for 
crystalline dust grains. Since the resolution of the IRAS images was approximately 2', the SCUBA 
data were also smoothed to this resolution. 

The fits indicate a fairly uniform dust temperature of ^ 21 it 2 K. There are no signs of 
strong temperature variations correlated with individual clouds, and this is further evidence that 
the brightness structure in the maps is dominated by column density variations. Therefore, we 
adopt a typical temperature of 20 K throughout the rest of this paper. 

If (3 is allowed to vary, the fit is better but the derived temperature docs not significantly 
change. Best-fit values of /? are /9 ~ 2.4, higher than the canonical value of /3 = 2. There is now 
much evidence for j3 > 2: Lis & Menten (1998) and Dowell et al. (1999) measured /3 ~ 2.8 and 
P ^ 2.5 respectively in GCMO. 25-1-0. 01. In other regions. Lis et al. (1998) measure /? values as 
high as 2.5 in the Orion Ridge and Bernard et al. (1999) also report PRONAOS observations of 
high (3 in the Polaris cirrus cloud. In laboratory work Koike et al. (1995) report (3 up to 2.7 in 
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carbon-based compounds at room temperature, and Agladze et al. (1996) have shown /3 > 2 in 
sihcates between 5 and 20 K. 

3.4. Spectral index map 

We have used our 450 /im and 850 fim maps to generate a map of spectral index a (defined by 
Su oc i^"), which is shown in figure 3. We first smoothed the 450 /um map to the lower resolution 
of the 850 nm observations, and imposed a lower threshold on both maps of ~ 3a, since spectral 
indices derived from small flux values are subject to large errors. 

Over most of the bright parts of the map a is in the range 3-4, indicating rather constant dust 
properties with a dust opacity index /3 ~ 2. The bulk spectral index for M-0.02-0.07 (the +50km/s 
cloud) is a = 3.5 ±0.05, whereas for M-0.13-0.08 (the +20km/s cloud) it is a = 3.7 ±0.05 (quoted 
errors are from the variation within the clouds). Assuming greybody emission at a dust temperature 
of 20 K. these indicate dust opacity indices of /? ~ 2.2-2.4. Although = 2 is the theoretical upper 
limit for crystalline dust grains, in section 3.3 we discuss some evidence for higher values of (3. 

Differences in spectral index may be due to differences in temperature. Since at these temper- 
atures, emission at 450 /v,m and 850 ;um is not truly in the Raylcigh- Jeans part of the spectrum, the 
difference between the spectral index a and the dust opacity index /3 is temperature dependent. 
It is therefore possible that two clouds may have the same /? but different temperatures, leading 
to different spectral indices. If this were the case, a lower a would indicate a cooler cloud. This 
cannot account for all the a variations: for example, a is low at only 3.2 ±0.1 for CO 0.02-0.02, 
and this cloud is thought to be relatively hot at ~ 60K (Oka et al. 1999). In section 3.3, we detect 
no significant temperature variations on ^ 2' scales. 

3.5. Cloud masses 

To estimate the total cloud masses from their integrated fluxes, we require the grain absorption 
cross-section per unit mass of gas and dust. We have adopted the parameters suggested by Pollack 
et al. (1994) with opacity index /? = 2, and a correction for Galactic Centre metallicity Z/Zq = 2 
(Mezger et al. 1989). We thus derive opacities of «;450/im = 4.94 x 10~^m^kg~^ and KssO/um = 
1.38 X lO^^m^kg"^. Assuming a uniform dust temperature of 20 K we have for A = 450 and 
850 //m, conversion factors of 63.3 and 513 Mq/Jj respectively. The total mass found by the 
survey, measured by integrating the total 850 /xm brightness, is then (53 ± 10) x IO^Mq. 

The largest source of error in these mass measurements is our uncertainty of the dust tempera- 
ture. The value of 20 K is derived in section 3.3, but Lis, Carlstrom, & Keene (1991) adopt a mean 
temperature of 30 K. However, the ISO study of GCM0.25±0.01 by Lis k Menten (1998) gives a 
mean temperature of ~ 18 K for this core, suggesting that the lower dust temperatures are correct. 
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We neglect optical depth effects as the emission is optically thin at 850 /xm except from Sgr B2. 

4. Further discussion 

The improved depth, resolution, and extent of the SCUBA survey reveal many new features in 
the Galactic Centre. Not only are there the known large star-forming cloud complexes but also a 
wide-ranging network of dusty filaments is revealed. For example, at negative Galactic longitudes 
close to 6 = there is filamentary structure which may be continuous and stretches along the 
Galactic plane for up to 0.8° (120 pc). There are partial shells, possibly associated with the known 
SNRs Sgr A East (GO.0+0.0), GO.3+0.0, GO.9+0.1, and Gl.4-0.1 (Green 1998). There is also a 
circular shell of dust clouds centred on l,b = (0.8, —0.18) which does not appear to be documented, 
and which we have labelled PPR GO. 80-0. 18. This may trace a wind-blown region or perhaps a 
supernova remnant. In a future paper, we will contrast these structures with star-forming regions 
near the Solar circle, comparing masses, densities, dust grain properties and star-forming efficiency. 

The total mass detected by our survey is (53 it 10) x 10® Mq, in the inner 400 pc. Dahmen 
et al. (1998) derived a weighted best estimate of the total mass in the inner 600 pc of the Galactic 
centre using a combination of dust emission, gamma rays, and C^^O molecular line emission. This 
estimate was 20-50 x 10® M©. The best optically-thin tracer in this set is the C-'^^O survey, from 
which was derived a mass of 17 x 10® in the central 600 pc. Although our survey covers a smaller 
area, the total mass detected is significantly greater. Therefore, the SCUBA data comprise the first 
optically-thin map to trace essentially all the mass in the CMZ at high resolution. 
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Fig. 1. — SCUBA Galactic Centre Survey: 450 /im flux density. 
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Fig. 2. — SCUBA Galactic Centre Survey: 850 /im flux density. 
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Spectral index map for central region, with contours of 850 nm map. 



